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ABSTRACT 


A general solution of the frequency equation for rigid rotors in undamped bearings is 
derived for forward and backward precession. The solution is applied to a wide range of 
speeds and shapes, with variation in center-of-gravity location, to obtain maps of the fre- 
quency characteristics and locate major and nonsynchronous critical speeds. When the 
center of gravity is at the bearing centerline midpoint, a set of two solutions exists. An 
additional set of two solutions appears, however, when the center of gravity is away from 
this point. No forward-precession high-frequency critical speeds exist if the moment- 
of-inertia ratio equals or exceeds the critical -speed ratio. With the center of gravity 
midway between single-row ball bearings and the polar and diametral moments of inertia 
equal, a nonsynchronous critical speed exists at all rotor speeds. 
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CRITICAL-SPEED ANALYSIS OF FLEXIBLY MOUNTED RIGID ROTORS 


by Richard H. Cavicchi 
Lewis Research Center 

SUMMARY 

A theoretical analysis of rigid rotors in undamped flexible bearings develops and 
solves a general frequency equation to determine frequencies of both forward and back- 
ward precession. Ibis study applies the general solution to a wide range of rotor rota- 
tional speeds and to shapes that vary from a pencil shape to a disk. The center of gravity 
varies from midway between the bearings to an outboard quarter point. Maps present the 
frequency characteristics of all these configurations for use as a guide in preliminary 
design. Besides locating major critical speeds, the maps locate numerous nonsynchro- 
nous critical speeds that may result from bearing defects. 

A single two-branch curve completely describes a set of two solutions for all shapes 
when the center of gravity lies midway between the bearings. With this location between 
single- row ball bearings and the polar and diametral moments of inertia equal, a nonsyn- 
chronous critical speed exists at all rotor speeds. When the center of gravity is away 
from the bearing centerline midpoint, two sets of solutions for forward and backward 
precession result. One is a low-frequency set, and the other is high. The frequency 
magnitude of forward precession exceeds that of backward precession in each set over 
all rotational speeds studied. As rotor speed increases, the frequency magnitude of for- 
ward precession increases and that of backward precession decreases. 

As the center of gravity is moved progressively from the bearing centerline midpoint, 
major critical speeds of the low-frequency set decrease in magnitude and those of the 
high-frequency set increase. The effect of the center-of-gravity location on major criti- 
cal speed decreases as geometries become more disk-like. 

This analysis revealed that, in backward precession, disks have lower major critical 
speeds than pencil shapes. No forward-precession high-frequency critical speeds exist 
if the ratio of polar to diametral moments of inertia equals or exceeds the critical- speed 
ratio. 



INTRODUCTION 


Numerous investigators have studied the problem of vibrations in rotating machinery. 
Gunter (ref. 1) traces the history of some of this work over the past century. He adds his 
own contribution with a general study of nonsynchronous precession in a rotor system. In 
references 2 and 3, Yamamoto makes major advances in this field by combining theoreti- 
cal and experimental treatments. He studies shafts and disks and varies the axial location 
and thickness of the disk. Yamamoto's work reveals that ball-bearing defects can cause 
critical speeds both below and above the major critical speed, the speed at which rotation- 
al and precession frequencies are equal. Dimentberg (ref. 4) also presents a variety of 
theoretical studies on shaft-disk systems backed by experiment. Eshleman and Eubanks 
report on a continuous shaft-disk system in reference 5. 

Despite the thorough treatment in these references, NASA requirements have stimu- 
lated further work on rotor systems. The work of reference 6 is aimed at predicting the 
deflection amplitude and frequency of vibrations in an axial turbopump, the design of which 
is described in this reference. 

The intention of the present report is to extend the prediction of vibration frequencies 
beyond the study of one particular machine as done in reference 6. By contrast, the work 
performed herein covers a wide range of rotor shapes and speeds. This report presents 
the results of a parametric analysis of rotor dynamic frequencies. Its purpose is to aid in 
preliminary design of turbomachinery by providing guidance as to the location of potential 
major and nonsynchronous critical speeds. 

In this study, the rotor geometry varies from a pencil shape to a disk, and the rotor 
speed is varied over a wide range. The analysis is performed for four axial locations for 
the center of gravity: midway between the bearings, at inboard and outboard quarter 
points, and at a bearing. . 

The results of calculations made for this analysis are presented in the form of dimen- 
sionless working charts relating vibration frequency with speed. Shape factors are used 
as parameters to display the effect of geometry on critical speed. This analysis is a con- 
tinuation of work suggested by Timoshenko in reference 7. 

The operating conditions of the axial turbopump investigated in reference 6 are located 
on some of the working charts to display their proximities to various critical speeds. 


DESCRIPTION OF ROTOR-BEARING SYSTEM 

Figure 1 is a representation of a rigid rotor on a shaft mounted in flexible bearings. 
The mass of the shaft is neglected. This study considers motion in the lateral (y and z) 
directions, but no axial motion (x). (All symbols are defined in appendix A. ) In figure 1, 
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y 



y- and z- displacements of the shaft are greatly magnified for clarity. This study further 
assumes no friction, no unbalance, and no external disturbing force. Specifically, there 
is no damping in the bearings. 

The rotor configuration varies from a pencil shape to a disk. Shape is introduced 
into the mathematical analysis in terms of the rotor polar moment of inertia I and its 
diametral moment of inertia 1^. To make the analysis nondimensional, the moments of 
inertia are expressed as dimensionless parameters: 


n 2 ~ 



( 1 ) 


w 3 



( 2 ) 


The values of these parameters are small for pencil shapes, and are larger for disks. 
For example, the ratio I/l^ for a rod whose length is 10 times its diameter is 0. 0149. 
For a disk of negligible thickness, it is 2. The parameter is the disk effect used in 
references 5 and 8. It is zero for a concentrated mass and is infinite for a disk having 
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all its mass distributed over a large radius. Therefore, geometries become more disk- 
like as 7 r g and increase. 

In the system studied in this report, the rotor center of gravity remains on the shaft 
centerline. The study does, however, investigate the effect of varying the axial location 
of the center of gravity. The disk effect v 2 and the moment-of-inertia ratio ir 9 each 
vary from zero to infinity in the maps prepared from this analysis. 


GENERAL ANALYSIS 


Figure 1 shows the geometry used in this study. The bearings are assumed to be 
linear springs, with a spring constant of k/2 for each bearing in both lateral directions. 


Equations of Motion 

Timoshenko presents the equations of motion for this analysis in reference 7 (pp. 287- 

288): 
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M 
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'?2 - yj) 
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W** 1 ♦**-«> 
1 a 2 2 

(3) 

kz 1 kz 9 

l.z 9 ) + A + z -0 
12 2 2 

(4) 

Z2-ZA _kl 2 z 2 ^ kl 1 z 1 _ 
If 2 2 

(5) 


Iu> 


^2 - *l \ , , / ?2 - , m 2?2 M l y l_, 


( 6 ) 


This study does not consider axial motion. Equations (3) and (4) are force equations; 
equations (5) and (6) are summations of torque. Equations (3) to (6) completely describe 
the free vibrations of a rigid rotor in flexible bearings. Timoshenko suggests the solu- 
tions 


yj = A sin pt 


(7) 
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y 2 = B sin pt 

(8) 

Zj = C cos pt 

0) 

z 2 = D cos pt 

(10) 


In equations (7) to (10), p is the frequency of rotor vibration. Calculation of p is the 
objective of this analysis. This study is not concerned with determining the deflections 
y and z, or the amplitudes A, B, C, and D. 

Substitution of equations (7) to (10) into equations (3) to (6) yields 



In matrix form, this set becomes 
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Frequency Equation 


Nontrivial solutions of A, B, C, and D in equations (11) to (14) require that the 
determinant of the coefficient matrix vanish. Expansion of the determinant and collection 
of similar terms yield 


(Ml/p 8 - 


M 2 I 2 w 2 + 2kMI 2 + kM 2 !^ 2 + ifj 


P 


2kMI 2 c 2 


/,2 , 2\ 2 3k MI 1 (,2 , 2 \ , 2 ,„ , , . . 2 t 2 

(^1 +l 2J — 2 vfl + V + k k *1 


, 2 t 2 2 k 3 M£ 2 

k I w + 


( 4 + 4 ) 


kVif 


9 i 4,4 
p 2 + M_ = 0 
16 


(16) 


Equation (16) is the general frequency equation for the model studied in this report. 


Nondimensional Frequency Equation 

The following dimensionless parameters make equation (16) nondimensional and 
shorten its appearance: 





7I 2 = 



(17) 


( 1 ) 


^3 



( 2 ) 


n 4 = 



(18) 
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( 19 ) 


<5 _ 2 

O = 7T j 


Substitution of these parameters into equation (16) yields 


2 C 4 
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2 2 2 
*2*3*4 





( 20 ) 


Solution of Frequency Equation 


Solution of the quartic equation (20) for S is inconvenient. It is simpler to solve for 
the shape-speed parameter The procedure is first to form a quadratic in i r 2 . 

Equation (20) takes the form 


2 [7,4 00 3 c 2 2 2 / 21 

7r g I S — 2S + S — ^ 3 ^ 4 ^ (S — 1) 1 





+ 


S 

2 



Factoring the coefficients of this expression yields 


S(S - l) 2 (s - 7 r§ 7 r|) 



( 22 ) 


It is convenient to solve this quadratic in for l/n^. 


The solution is 
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(23) 


4(S - 1) 




5 ^* 2^4 



The result of solving equation (23) for yields the general solution of the frequency 

equation 

2 


4^S (S - 1)jt 2 


6^371^ 




!i\ + if 


s - 1 


(24) 


In equations (23) and (24), for forward precession, 

6 = 1 


and for backward precession, 


6 = -1 


(25) 


(26) 


The presence of the factor (S - 1) in the denominator reveals that no solution exists 
for 


S = 1 


or 



8 


Major Critical Speed 

A major critical speed is said to occur when the magnitudes of the precession and 



rotational frequencies are equal. Also called synchronous speed, this condition is repre- 
sented mathematically by 


n l 

r cr = 7 = ±1 


or 


(27) 


P =±1 

co 


When this relation is used in the general solution of the frequency equation (eq. (24)), 
the major critical speed results as follows: 


-I 1/2 
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2 
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(28) 


For backward precession (6 = -1), a positive sign before the radical in equation (28) al- 
ways yields the high-frequency solutions; a minus sign yields the low-frequency set. For 
forward precession (6 = 1), the same pattern is true when < 1. When tt ^ > 1, how- 
ever, the effects of the signs are reversed. 


Calculation Procedure 

Ihe calculation procedure uses S as an independent variable for which a range of 
values is assigned. This is an indirect method that evaluates the dimensionless frequency 
parameter from graphs. From equations (17) and (19), 
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77 , = 


P 






Appendix B applies equation (24) to the four center -of- gravity axial locations. The 
equations for these four cases are collected here for Comparison. 


Center of gravity midway between bearings, case 1 ( l^/l = 1/2): 


6 73 " , 


1*4 = — F~ 

4 Vs vr 2 


(B2) 


Center of gravity at quarter point, case 2 = 3/4): 


5 V4 = V® 


5S _ x 


4^S (S - 1)7 Tc 


(B15) 


Center of gravity at one bearing, case 3 = 1): 


»v 4 =V®- 


2S - 1 


4^S (S - 1)7 


(B17) 


Center of gravity at outboard quarter point, case 4 - 5/4): 


13S 


577^77^ = ^/S - 


- 1 


4^S (S - 1)77, 


(B19) 


The effect of the axial location of the center of gravity appears in these four equations and 
in equation (24) only as a coefficient of S, 2 • Tilis factor is designated the 

center- of- gravity influence coefficient. Thus, equation (24) is a relatively weak function 
of the center- of- gravity location. The center- of- gravity influence coefficient varies from 
4/4 to 13/4 as the center- of- gravity location varies from midway between the bearings to 
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an outboard quarter point l^/l = 5/4. Figure 2 shows the variation of this coefficient 
with the center- of- gravity axial location. 

Equation (24) can easily be applied to cases other than the four presented in this re- 
port. For quick estimates, however, the graphs of equation (B2) can be used for center- 
of-gravity axial locations with 1/2 ^ |Zj/j| < 3/4, since figure 2 reveals that the in- 
fluence coefficient does not change greatly in this region. 

Equations (B15), (B17), and (B19) reveal that there is no solution for S = 1 for 
cases 2, 3, and 4. For case 1, however, all solutions exist. 



Center-of-gravity axial location, 


Figure 2. - Variation in center-of-gravity influence coefficient with 
center-of-gravity axial location. 
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DISCUSSION AND RESULTS 
Universal Plots 


The basic graphs of this report are the universal plots. They apply to a wide range 
of geometry and speed. The universal plots are made from equations (B2), (B15), (B17), 
and (B19). 

Appendix B shows that a single two-branch curve can represent all geometries and 
speeds for case 1 in which the center of gravity is midway between the bearings. This 
abbreviated universal plot (made from eq. (B4)), shown in figure 3, has the adjusted 
shape-speed parameter 




OJ 



< 



Figure 3. - Abbreviated universal plot for center of gravity midway between 
bearings, case 1; equation (B4). 
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as the abscissa and the adjusted-frequency parameter 


P 



as the ordinate. The upper branch in figure 3 represents forward precession, and the 
lower branch, backward precession. When 6 = 1, positive values of ^3^4 in equa- 
tion (B4) result from positive values of 7 to yield the upper curve in figure 3. When 
6 = -1, positive values of ^3^4 result from negative values of 77^ to yield the lower 
curve. In reference 8, Den Hartog observed, presumably from experiment, that the low- 
frequency branch represents backward precession. 

The backward- precession branch should be shown in the fourth quadrant, since it rep- 
resents a direction opposite from the upper branch. Both branches are shown in the first 
quadrant to save space, however. Figure 3 displays only one-half of the total solution of 
equation (B4). Although not shown, a mirror image of figure 3 in the second quadrant 
represents the other symmetrical half. 

Because equation (16) is of the eighth order in p, there are four additional solutions 
for case 1. These four solutions are represented by 

(k - Mp 2 )(k - Mp 2 ) = 0 (29) 

o 

This equation results from setting l^/l = 1/2 which permits factoring (k - Mp ) from the 
first two rows of the matrix in equation (15). Each factor in equation (29) yields 

= ±1 (30) 


Yamamoto states in reference 2 that the precession frequency is constant when the 
center of gravity is midway between the bearings. His constant value is represented by 
equation (30), which, as discussed previously, represents only the less-interesting half of 
the eight solutions. Yamamoto's presentation in reference 2 therefore neglects the vari- 
able solutions of figure 3 and equation (B2) of this report. 

For the axial turbopump studied in reference 6, the product 
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is 0. 19 at design speed. This value for the abscissa in figure 3 locates the reference 
turbopump design operating point shown by the open circles. 

When the center of gravity does not lie midway between the bearings, the universal 
plots consist of a family of curves with 7 t 2 as the parameter. Figure 4 presents univer- 



Shape-speed parameter, 

(a) Center of gravity midway between bearings-, case 1; equation (B2). 
Figure 4. - Universal plots. 









Shape-speed parameter, 7^4 = (i/I^&yVS/M) 

(d) Center of gravity at outboard quarter point; case 4; equation <B19), 
Figure 4. - Concluded. 





sal plots of this type for all four cases. In these figures the shape-speed parameter 


I <jJ 



is the abscissa, and the ordinate is the frequency parameter 


P 



Figure 4(a) repeats the information for case 1, presented more concisely in figure 3, to 
provide a direct comparison of case 1 with the other three cases, as represented in fig- 
ures 4(b) to (d). It is not possible to reduce the family of curves to one two-branch curve 
for cases 2 to 4. It may be of interest to note that the set of curves in figure 3 corre- 
sponds to the set in figure 4(a) (not shown) for = 0. 25. This fact results from equating 
the corresponding ordinate or abscissa parameters (2 V ^2 *1 = ff l>* 

Figures 4(b) to (d) each present two forward-precession solutions (solid curves) and 
two backward (dashed curves). These four solutions plus their mirror images in the 
second quadrant (not shown in fig. 4) make up the eight solutions called for by equa- 
tion (16). 

Intersections of the curves with the ordinate axis of the universal plots represent vi- 
bration frequencies of nonrotating rotors. There is no distinction between forward and 
backward precession at these points, as the two branches of the curves join there. The 
solution given by equation (30) is of a different mode from the zero- speed solutions given 
in figure 4(a) by the ordinate-axis intersections. 

The universal plots in figures 4(b) to (d) display the no-solution condition of equa- 
tion (24) when p/^k/M is plus or minus 1. A no-solution barrier does not exist, of 
course, for case 1, as shown in figure 4(a). For cases 2, 3, and 4, the no-solution bar- 
riers divide the plots into a low-frequency set and a high-frequency set. There is a 
forward- and a backward-precession branch in each set. 

As rotor shapes become more disk-like and the disk effect tt 2 approaches infinity, 
curves on the universal plots approach a 45° line through the origin for all four cases. 

This condition is obvious from equation (23). It occurs both below and above the no- 
solution barrier of figures 4(b) to (d). In fact, equation (23) also shows that the no- 
solution barrier itself (tt^ = 1) also represents geometries in which is infinite. 

In figures 4(b) to (d), forward precession curves in the low-frequency set all cut the 
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45° line through the origin in a common point. This observation can be shown analytically 
by setting 

6 = 1 


and 


*3*4 = V® = n l 


in equation (24). The result is 



( 31 ) 


which is independent of tt 2 . 

Comparison of the universal plots in figure 4 reveals that, for a given shape and ro- 
tational speed, the effect of the center-of-gravity axial location as it is displaced outward 
from the midway point is a small rise in the magnitude of precession frequency. 


Frequency Plots 

Curves of precession frequency plotted against rotational frequency are helpful in 
visualizing the locations of critical speed for a given geometry. Figure 5 presents a few 
frequency plots for selected geometries of cases 1 and 3. Less general than the universal 
plots, the frequency plots are obtainable from them or from equations (B4) and (24) by 
specifying tr^. When is 1, frequency plots and universal plots are identical. 

Each curve in figure 5 denoting constant represents a constant geometry. The 
design operating point of the reference turbopump is shown by the two lowest open circles 
in figure 5(a). The straight lines in figure 5 denote various critical speeds, to be discus- 
sed in the following two sections. Thus, each curve in figure 5 traces the passage of a 
given machine through potential critical speeds as it accelerates or decelerates. 
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(a) Center of gravity midway between bearings; case 1. 
Figure 5. - Frequency plots showing critical speeds. 






Frequency parameter, 7Tj ■ [/v U/M 







Frequency parameter, 







Major Critical Speeds 


An important aspect of this investigation is to locate the major critical speeds for all 
geometries considered. This information is significant because forward-precession ma- 
jor critical speeds may be excited by rotor unbalance. For this purpose, plots relating 
major critical speed with geometry are helpful. Major critical speed plots, presented in 
figures 6 and 7, are derived from equation (28). They may also be obtained from inter- 
sections of the curves on the universal plots with straight lines through the origin having 
slopes of l/Vg = l^/l. 

The abbreviated major critical- speed plot for case 1 consists of one two-branch 
curve and is shown in figure 6. The abscissa in this figure is lAp and the ordinate is 



The curves in figure 6 are derived from equation (B9). The upper branch of the curve 
gives the solution for forward precession; the lower branch is for backward precession. 



Figure 6. - Abbreviated major critical-speed plot for center of gravity 
midway between bearings, case 1; equation (B9). 
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Figure 7. - Major critical-speed plots. 


Thus, for case 1, major critical speeds, when they occur, are higher for forward than 
backward precession for any shape. Equation (B9) and figure 6 show that, for constant 
7Tg for case 1, the major critical speed is inversely proportional to yir^. 

Equation (B9) shows that, for forward precession (6 = 1), the major critical speed 
becomes infinite when is 1 and is imaginary beyond 1. The upper curve in figure 6 
illustrates these conditions. This fact implies that forward major critical speed can be 
avoided entirely for case 1 by using designs having I/l^ greater than 1. If such designs 
are not feasible, designing with I/I^ as close to 1 as possible may yield a major critical 
speed well above the design speed. 

Because the dashed curve in figure 6 remains finite, backward-precession major 
critical speeds cannot be so easily avoided by selective design. However, if designs with 
a low value of ir 2 (pencil shapes) are used, the major critical speed may be high enough 
to exceed design speed. 

The steep slope of the upper curve in figure 6 reveals that the major critical speed of 
forward precession is very sensitive to shape. The dashed curve shows that the major 
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Figure 7. - Continued. 


critical speed of backward precession is much less sensitive. Also, the adjusted major 
critical-speed parameter for backward precession decreases with increasing 1/1^. 

According to the locations of the lowest two open circles in figure 5(a), the axial 
turbopump of reference 6 passes through a major critical speed as it accelerates to the 
design point. Major critical-speed plots for all four cases are shown in figure 7. These 
figures have Ij/M2 2 as the abscissa and w cr / Vk/M as the ordinate, with I/E^ as a pa- 
rameter. Figure 7(a), which repeats the information given in figure 6, aids in a compari- 
son. As in the universal plots, there are two families of curves on the major critical- 
speed plots of cases 2 to 4. One family is above the no-solution condition, and one is be- 
low. As noted earlier, the positive sign before the radical in equation (28) yields the high- 
frequency set. The curves in figure 7 show that, in general, major critical speed de- 
creases as Ij/MZ^ increases. For cases 2 to 4, figures 7(b) to (d) show an exception to 
this observation for low-frequency forward precession when I/l^ is greater than l. When 
I/Ij is 1, the forward major critical speed is independent of I^/MZ 2 for cases 2 to 4. 

This observation arises from setting 6 = 1, n^ = ir^, and 77 ^ = 1 in equation (23). The 
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(c) Center of gravity at one bearing; case 3. 
Figure 7. - Continued. 


result is again equation (30), which was derived earlier for a different situation. 

On the other hand, for forward precession, major critical speed increases as I/Ij 
increases. The opposite variation occurs for backward precession. Thus, for backward 
precession, disks have lower major critical speeds than pencil-shape configurations. No 
such conclusive statement can be made for forward precession because of the opposite ef- 
fects of Ij/M l 2 and l/l v 

At values of I. /Ml 2 greater than about 1, its effect on major critical speed is neg- 
ligible for the high-frequency set. Furthermore, in this range of Ij/MZ , variations in 
1/1^ have a minor effect as the curves for this set converge. 

It was observed for case 1 from equation (B9) and in figure 6 that no major critical 
speed exists for forward precession when I/lj is 1 or higher. This observation is true 
for the other three cases, too, but only for forward precession in the high-frequency set. 
Equation (28) is useful in explaining this circumstance. When 6 = 1 and ir^ = 1, it is 
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Figure 7. - Concluded. 


evident from equation (28) that co cr / yk/M becomes infinite. When = 1 and 7 Tg > 1, 
the magnitude of the expression following the ± sign exceeds that preceding it. There- 
fore, c cr /yk/M from equation (28) is imaginary when the minus sign is used; this con- 
dition pertains to the high-frequency set. 

Consideration of the universal plots may give added insight to this explanation. A 45° 
line through the origin on a universal plot not only locates the major critical speeds for an 

O 

I/lj of 1, but also represents a value of infinity for I^/M l . This condition is also evi- 
dent from equation (24) when is infinite. The forward-precession curves of the high- 
frequency set therefore cannot intersect the 45° line when I/I^ is 1. Thus, no major 
critical speeds are possible for forward precession of the high-frequency set when l/l^ 
is 1 or greater. Figures 3 and 4(b) to (d) all confirm this observation. It is apparent, 
however, that intersections with the forward-precession low-frequency curves still occur. 

Figure 8 presents the major critical-speed results in a form to show the effect of 
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Precession 

Forward 

Backward 



Figure 8. - Continued. 




Figure 8. - Continued. 


axial location of the center of gravity. Four plots are presented in figure 8, each for a 

9 

selected value of l/l j. In addition, I^/M l varies as a parameter on each plot. Hie 
plots in figure 8 are derived directly from the corresponding plots of figure 7 or from 
equation (28). 

To avoid confusion among the many curves, the backward-precession curves in fig- 
ure 8 are shown in the fourth quadrant. As usual, they are dashed. Figure 8 has a mir- 
ror image in the second and third quadrants. The abscissa range of l^/l of the image 
is from -0. 25 to 0. 5; that is, ^/l varies from 0. 5 to 1. 25. 

There are two families of curves for both forward and backward precession except 
for l^/l of 0. 5, as noted in the discussion of figure 4. The merging curves in figure 8 
graphically illustrate the fact discussed in connection with figure 3 that the two sets of 
solutions ostensibly absent for case 1 actually exist in = 1. The slopes of all the cur- 
ves in figure 8 are zero at the bearing centerline midpoint because of symmetry. Differ- 
entiating equation (28) and using l */Z = 1/2 proves this observation mathematically. 
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Figure 8. - Concluded. 


Figure 8 shows that the effect of the axial location of the center of gravity decreases 
as both 77^2 and increase. Therefore, the effect of the axial center-of-gravity loca- 
tion decreases as geometries become more disk-like. 

Furthermore, there are no forward precession high-frequency curves in figures 8(c) 
and (d) for reasons mentioned on pages 27 and 28. 

In figure 8, solutions that lie within the - 1 to 1 ordinate band reveal a trend of de- 
creasing major critical-speed magnitude as the center of gravity is located progressively 
outward from the bearing centerline midpoint. Both families outside this band display the 
opposite trend. 

In figure 8(c), in which 7!^ = 1, the forward-precession low-frequency curves are 
identical for all values of This characteristic, observed in figures 7(b) to (d), re- 

sults from equation (31). 

Figures 7 and 8 should aid in preliminary design by pointing out potential major cri- 
tical speeds for various configurations. 
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Nonsynchronous Critical Speeds 

It is well known that critical speeds other than the major critical speed may induce 
destructive vibrations. The accompanying motions of the rotor are commonly called non- 
synchronous precessions. 

Yamamoto has made the most thorough investigation of nonsynchronous precessions 
in rotors supported by ball bearings. He has determined that the bearings are respon- 
sible for this phenomenon. In reference 2, he reports that two classes of nonsynchronous 
critical speeds occur when double- row self-alining ball bearings are used. These mo- 
tions occur at rotor speeds both greater and less than synchronous critical. The spring 
characteristics are linear for this bearing type. Therefore, the rotors studied herein 
may be assumed to be mounted in double-row ball bearings because of the linear spring 
constant used in this analysis. 

Yamamoto's work in reference 3 reveals two other classes of nonsynchronous criti- 
cal speed for single-row ball-bearing configurations. These motions all occur at rotor 
rotational speeds above synchronous. This type of bearing has nonlinear and nonsymmet- 
rical spring characteristics because of the absence of self alinement in the bearings. 

Therefore, the conditions for nonsynchronous critical speeds resulting from single- 
row ball-bearing defects cannot be rigorously applied to the solutions of the present anal- 
ysis. However, even for double-row ball-bearing defects, the present analysis is approx- 
imate because of its idealized assumption of no damping. Furthermore, Yamamoto shows 
in reference 3 that theoretical solutions from his linear spring- constant analysis are good 
approximations to his experimental results obtained with single-row ball bearings (non- 
linear spring constants). Therefore, nonlinear spring constant conditions are applied 
herein without hesitation to the present analysis for the purpose of locating potential non- 
synchronous critical speeds. 

Double-ro w ball bearings . - In experiments with self-alining ball bearings, 

Yamamoto observed that very slight defects can produce driving forces at frequencies dif- 
ferent from the shaft speed. One class of such nonsynchronous motion arises from non- 
uniformity in the ball diameters. Another set of nonsynchronous precession may result 
from noncircular inner or outer bearing races. Yamamoto observed that all commercial 
bearings made in Japan, and presumably elsewhere, have these defects. 

Nonuniform ball diameters cause nonsynchronous precession at rotor speeds higher 
than critical. Yamamoto's work disclosed that the most serious conditions of this type 
occur when the critical- speed ratio is 


r =£. = 

u> 2. 65 


(BIO) 
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for forward precession. The constant in this expression is an average value obtained 
experimentally and theoretically by Yamamoto, and reported in references 2 and 9. He 
presents the following relation in these references: 

r = (32) 

2 + 2 i 

For backward precession, Yamamoto finds 

r = — = — — (Bll) 

w 4. 1 

Appendix B uses these expressions in equation (B8) for case 1 yielding equations (B12) 
and (B13). These equations are plotted in figure 9, which is a companion plot of figure 6, 
with appropriate change of constant. Figure 9 shows that forward nonsynchronous pre- 
cession of this type cannot occur for shapes having I/l^ greater than the critical- speed 
ratio 1/2. 65. Thus, forward-precession nonsynchronous speeds of this type are limited 
to a small range of ir^. Backward precession is possible over the entire range of I/I^, 
however. 

At this point, the use of equation (B8) is convenient in making a general statement: 

2 V^2 77 A = - -- - (B8) 

Y r(r-5 77g) 

For case 1, no forward-precession critical speeds (either synchronous or nonsynchronous) 
exist when the moment-of-inertia ratio ir^ equals or exceeds the critical-speed ratio r. 

Yamamoto's work revealed that the amplitudes during backward-nonsynchronous mo- 
tion are smaller than in forward precession. Therefore, although the dashed curve re- 
mains finite in figure 9, the small amplitudes reduce the severity of the consequences of 
this motion. 

Lines representing r = 1/2. 65 and r = 1/4. 1 for case 1 are shown in the frequency 
plot of figure 5(a). Likewise, corresponding lines are shown in the frequency plots of 
figures 5(b) to (d) for case 3. 

In figure 5, only the intersections of the r = 1/2. 65 lines with the forward- 
precession curves are meaningful. Similarly, the meaningful intersections for the 
r = 1/4. 1 lines are with the backward-precession curves. Solid circles in figure 5 de- 
note valid intersections. 
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Figure 9. - Abbreviated nonsynch ronous critical-speed plot for 
r = 1/2. 65 and r = -1/4. 1 for center of gravity midway 
between bearings; case 1, equations (B12) and (B13). 


Most of the r = 1/2. 65 and r = 1/4. 1 intersections in figures 5(b) to (d) occur 
below the no-solution barrier |p/y(k/M) | = 1 for both forward and backward preces- 
sion. Although these intersections are at the lower end of the rotor-rotational-speed 
scale, they occur at more than twice the major critical speeds. 

In reference 2, Yamamoto points out that noncircular inner and outer bearing races 
cause nonsynchronous forward precession at critical- speed ratios of 2, 3, and 5. The 
r = 5 motion is not serious, however, because of the low amplitudes he observed at this 
condition. Lines denoting critical- speed ratios of 2 and 3 are shown on the frequency 
plots of cases 1 and 3 in figure 5. All intersections, shown by solid circles, take place 
at rotor rotational speeds well below the major critical speeds. This observation im- 
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plies that even low-speed machines can experience this type of nonsynchronous preces- 
sion. Furthermore, most machines, including the reference turbopump, must acceler- 
ate through these critical speeds to reach design speed. 

Figure 5(d) shows that when l/l^ is 2, no high-frequency forward-precession inter- 
sections exist with the r = 2 lines. This circumstance is analagous to that for case 1 in 
which it was noted that, by equation (B8), no forward precession exists when u g > r. 

For cases 2 to 4, both a graphical and a mathematical explanation aid in understanding 
the absence of forward-precision intersections when 7ig = r = 2. 

Although not presented herein, frequency plots for cases 2 to 4 made for 1/1^ of 2 
would also disclose no nonsynchronous critical speeds for the forward-precession high- 
frequency sets when r = 2. Consideration of the location of the = 00 line on a fre- 
quency plot will illuminate this situation because no forward-precession high-frequency 
solutions are allowed below this line. For example, on a frequency plot constructed for 
7Tg = 2, the ^ ne c °i nc id es with a critical-speed-ratio line of r = 2 and therefore 

allows no high-frequency solutions. More generally, coincidence of the Wg = 00 line and 
a critical-speed- ratio line for forward precession (6 = 1) on a frequency plot for ir^ = r 
can be proved mathematically from equation (23). 

Therefore, a general conclusion can be stated for all cases. Forward-precession 
high-frequency critical speeds do not exist for geometries in which I/I^ equals or ex- 
ceeds the critical-speed ratio r. Figures 4 and 5 show that forward-precession low- 
frequency intersections exist, however. 

It is unlikely that a set of bearings would have either nonuniform ball diameters or 
noncircular bearing races without having both defects. Thus, it may be assumed that all 
the critical speeds revealed for both kinds of bearing defect will occur. 

Single- row ball bearings. - One class of nonsynchronous motion associated with 
single-row ball bearings is called subharmonic. Yamamoto observed this motion only in 
forward precession. Subharmonic motion is characterized by critical-speed ratios r of 
1/2, 1/3, . . . 1/n, where n is a positive, real integer. It can occur only after the 
major critical speed is exceeded. 

Lines for subharmonic motion are shown on the frequency plots for cases 1 and 3 in 
figure 5. Only intersections with forward-precession curves are valid, as shown by solid 
circles. Figure 5(a) shows that no high-frequency nonsynchronous precession of this type 
is possible for case 1 when 1/1^ is 0. 5 or greater. This observation is also true for all 
cases, and is a consequence of the general conclusion stated earlier relating n g and r. 

The axial turbopump studied in reference 6 was mounted in single-row ball bearings. 
According to figure 5(a), the locations of the two lowest open circles reveal that it 
passes through the r = 1/2 nonsynchronous critical speed when accelerating to the 
design point. In fact, the design point itself lies close to this nonsynchronous condi- 
tion. 
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Yamamoto further observed that when single- row ball bearings are used, two natural 
frequencies build up together at certain speeds. The absolute value of their sum or dif- 
ference is related to rotor rotational speed by 

<u=|p i ±p.| (33) 

in which p^ and Pj are any two of the natural frequencies as calculated from equa- 
tion (24). In reference 3, Yamamoto reports that vibration amplitudes from these sources 
may exceed those of the major critical speed. It is evident from any of the universal 
plots that the condition i = j is not allowed, except for the trivial case of zero rotor ro- 
tational speed. In this instance, only 


Pl = P 4 


and 


P 2 = 03 


are allowed. 

It was observed in figure 5(a) that there are only two solutions for case 1 for each 
value of I fly The only summed-and-differential motion possible for this case is there- 
fore 

= |Pi " P4I 

where p^ represents forward-precession solutions, and p^ represents backward- 
precession solutions. 

The uppermost open circle in figure 5(a), obtained by adding the ordinates of the two 
lowest open circles, represents the summed-and-differential motion for the reference 
axial turbopump. The location of this circle relative to the 45° line through the origin 
shows that this turbopump accelerates through a summed-and-differential nonsynchro- 
nous critical speed just prior to reaching its design point. 

Figure 10(a) presents a copy of the = 0. 4 curve of figure 5(a) with curves of 







Frequency parameter, 7^ = p/VkTM 



Rotational -speed parameter, 7 r 4 = uj/Vk/M 


<b) Center of gravity at one bearing; case 3; 7r 2 = * 0. 1; 7r 3 = I/Ij = 0.4. 

Figure 10. - Concluded. 
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added to illustrate the possibilities for case 1. The uppermost solid circle locates a non- 
synchronous critical speed at 2. 16 for 



This point is established by the intersection of the 


PfP 4 



curve with the 45° line through the origin. It is clear that no intersection is possible for 
the , , 


Pj +P4 



curve. 

Either equation (B4) or figure 5(a) can be used to show that |p^ + P4I coincides with 
the line | (pj ± Pj)|/o> = 1 when I/lj has a value of 1. From this observation, it is evi- 
dent that designs having the center of gravity midway between the bearings and I = Ij 
should be avoided if single-row ball bearings are used because all rotational speeds are 
critical. 

For a rotor with four natural frequencies, as observed for cases 2 to 4, there are 
twelve summed-and-differential combinations. Figure 10(b) presents a copy from fig- 
ure 5(c) of the four solutions for case 3 when 772 = 0. 1 and 773 = 0. 4. Figure 10(b) also 
displays all 12 possible |p. ± p^| combinations. 

Because at zero rotor-rotational speed p^ coincides with p^ and Pg coincides with 
Pg, no IP4+P4I or Ipg+Pgl intersections with the 45° line through the origin of a 
frequency plot are possible. Figure 10(b), which is drawn for shapes intermediate between 
pencil shapes and disks, shows 10 |p^±pj| intersections with the line r = 1. There- 
fore, the maximum number of intersections for any geometry is 10. 

Yamamoto and Hayashi have proved mathematically in reference 10 that only summed 
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harmonic oscillations can occur if both natural frequencies Pj and p^ represent forward 
precession or if both represent backward precession. If one frequency represents for- 
ward precession and the other, backward, only differential harmonic oscillations can oc- 
cur. The following tabulations summarize the combinations that can and cannot occur. 

Can occur: 

CO = 

c o = 

CO = 

CO = 

CO = 

CO = 

Cannot occur: 

CO = 

CO = 

CO = 

CO = 

CO = 

CO = 

Valid intersections with the major critical- speed line, according to the allowed combina- 
tions of equation (34), are shown by the solid circles in figure 10(b). The four major cri- 
tical speeds also are located by the appropriate solid circles. 

For disks, no intersections for jpj* Pj | can exist, as figure 5(d) will verify. There- 
fore, there are five |p. ± pj | intersections with the line r = 1 on a frequency plot for a 
disk. But equation (34) shows that only |p 2 “ P 3 I > 1 Pg ~ P4 1 » anc * 1 P3 + P4I can occur 


iPl + P 2 1 
IPl- P 3 | 

|Pl ‘ P4I 
|P2 " P3I 
|P 2 ‘ P4! 
IP3+P4I J 


(34) 


|Pl " P 2 I 
!pi + p 3 I 

IPl + P 4 I 
!p 2 +P3I 
lp 2 + P 4 I 
IPs - P4 1 


-\ 


(35) 
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for disks. In references 3 and 10, Yamamoto presents experimental evidence to substan- 
tiate equations (34) and (35) as applied to disks. 

A conclusion from this phase of the study is that the maximum number of summed- 
and-differential nonsynchronous critical speeds varies from three for a disk to six for 
pencil shapes. 

Plots similar to figure 8 may be drawn for any of the nonsynchronous critical speeds 
discussed, as well as for the major critical speeds. 


SUMMARY OF RESULTS 


This rotor critical- speed analysis has yielded the following results when the center of 
gravity is midway between the bearings: 

1. The solution for rotor frequencies simplifies to a single two-branch curve, one 
for forward and one for backward precession. 

2. The forward major critical speed rises rapidly and becomes infinite as the ratio 
of polar to diametral moments of inertia I/Tj approaches 1. No forward-precession ma- 
jor critical speeds exist when I/lj equals or exceeds 1. The magnitude of backward 
precession frequency, however, remains finite and decreases slowly throughout the range 

of l/l v i - 

3. Major critical speeds are inversely proportional to l for constant I/I^. 

4. When single- row ball bearings are used, designs should avoid having the polar and 
diametral moments of inertia equal because a summed-and-differential nonsynchronous 
critical speed occurs for all rotor rotational speeds. 

5. An experimental axial turbopump passes through a major and a nonsynchronous 
critical speed during acceleration, and at its design point operates close to nonsynchro- 
nous critical speeds. 

For more general center -of -gravity locations, the following observations were made: 

6. A low-frequency set of solutions occurs for forward and backward precession 
within a band width of - 1 to 1 for p/^k/M. A high-frequency set occurs outside this 
band. A no-solution barrier exists at ±1 for p/yk/M except when the center of gravity 
is centrally located. 

7. At any rotor speed of a given configuration, the frequency magnitude of forward 
precession exceeds that of backward precession. 

8. The frequency magnitude of forward precessions increases with increasing 
rotor speed. The opposite trend occurs for backward precession. 

9. As the center- of- gravity location is moved progressively outward, major critical 
speeds of the low-frequency families decrease in magnitude and those of the high- 
frequency families increase. 
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10. In backward precession, disks have lower major critical speeds than pencil 
shapes. 

11. No forward-precession high-frequency critical speeds exist for geometries in 
which the polar to diametral moment- of- inertia ratio equals or exceeds the critical-speed 
ratio. 

12. The effect of the axial location of the center of gravity on the major critical 
speeds decreases as geometries become more disk-like. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 6, 1968, 

122-29-02-20-22. 
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APPENDIX A 


SYMBOLS 


A, B, C, D 

vibration amplitude, m 

r 

critical-speed ratio, 

2 

outer diameter of inner 


^ 1 A 4 = p/w 


ring of ball bearing, m 

S 

square of frequency pa- 

d 

diameter of individual ball 
in ball bearing, m 


rameter, (p/yk/M)^ 

I 

rotor polar moment of 

t 

time, sec 


inertia, kg-m^ 

X 

axial coordinate of arbi- 

I 1 

rotor diametral moment 


trary reference frame, 

1 

2 

of inertia, kg-m 


coincident with shaft 
centerline when undis- 

k 

linear spring constant, 
kg/sec ^ 

y 

turbed, m 

lateral coordinate perpen- 

yk/M 

natural frequency used for 
nondimensional pur- 
poses, rad/sec 

yi,y 2 

dicular to x- coordinate, 
m 

lateral distance from ar- 

z 

distance between bear- 
ings, m 

1 Li 

bitrary reference frame 
to bearings in y- 

Zj, Zg 

distances from bearings to 


direction, m 

rotor center of gravity, 
m 

Z 

lateral coordinate mutually 
perpendicular to x- and 

M 

rotor mass, kg 


y- coordinates, m 

n 

general factor in har- 
monic sequence 
1, 1/2, 1/3, ... 1/n 

z l’ z 2 

lateral distance from ar- 
bitrary reference frame 
to bearings in z- 
direction, m 

P 

rotor precession fre- 



quency, rad/sec 

6 

±1 

P 1’ P 2’ P 3’ P 4 

rotor precession frequency 
used in summed-and- 
differential nonsynchro- 
nous critical speeds, 
rad/sec 


frequency parameter, 
p/^k/M = 
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*2 

dimensionless rotor moment of 

Subscripts: 


inertia or disk effect, I^/MZ^ 

B 

backward 

n 3 

ratio of rotor polar to diametral 
moment of inertia, I/I^ 

cr 

major critical 


F 

forward 

TTa 

rotor rotational speed param- 



TX 

eter, u>/ ^/k/M 

ij 

dummy indexes, 

CO 

rotor rotational speed, rad/sec 

ns 

nonsynchronous 


i - 1, 2, 3, 4 
j = 1, 2,3,4 
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APPENDIX B 


APPLICATIONS 

This appendix applies the general solution to the frequency equation (24) to four spe- 
cific cases. The only variation in these cases is the axial location of the center of grav- 
ity. 


Center of Gravity Midway Between Bearings (Case 1) 

The simplest case is that in which the center of gravity lies midway between the bear- 
ings. For this condition (see fig. 1) 



1 

2 


(Bl) 


Insertion of this relation into equation (24) gives 

87r 3 77 4 = ys L_ (B2) 

It is of great significance that the factor (S - 1) in equation (24) cancels for this case. 
Basically, it reduces equation (24) from second order in S to first order. By equa- 
tion (19), equation (B2) is therefore a quadratic in v j. Rewritten with ir^, equation (B2) 
is 


5w 3 7T 4 ~ n l 


1 _ 

4ff l*2 


(B3) 


If this expression is plotted with n ^ and as coordinates, a family of curves with 

^2 as the parameter results. More judicious choice of coordinates, however, reduces 
the family of curves to a single curve. The procedure is to use in each of the two 
coordinates. Thus 



(B4) 


46 



The coordinates used are 


and 



(B5) 


(B6) 


Den Hartog suggested this technique in reference 8 (p. 259). 
The critical speed ratio is defined by 


if* 

r = — (B7) 

n 4 


With this relation, equation (B4) becomes 


2 V^2 n 4 = ' — (B8) 

yr(r - Swg) 


For the major critical speed, defined by equation (27), r = r cr = 1, equation (B8) becomes 


2 V"2 l, 4,cr° 1 <B9) 

V 1 - 6 ’3 

Equation (28) also reduces to this form when equation (Bl) is used. 

In reference 2, Yamamoto observed the following nonsynchronous critical speeds, 
among others: 

Forward precession: 

r = — L_ (BIO) 

2. 65 
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Backward precession: 


r = -- 


4. 1 


Use of these expressions and the appropriate value of 6 in equation (B4) yields 
Forward precession: 



Backward precession: 



Center of Gravity at Quarter Point (Case 2) 

If the center of gravity is located one-fourth of the axial distance between the 
ings, figure 1 shows that 


L _1 = 3 
l 4 


l 


1=1 

l 4 


With these relations, equation (24) becomes 


5S 


8 *3*4 = Vs- 


4^S(S - l)vr 2 
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(Bll) 


(B12) 


(B13) 


bear- 


(B14) 


(B15) 


i in i i 


The change-of-coordinate technique used in equation (B4) is of no value in this case 
because of the S in the numerator of equation (B15); that is, there is no choice of coor- 
dinates that can reduce the family of curves to a single curve. 


Center of Gravity at One Bearing (Case 3) 


When the center of gravity is at one of the bearings, the conditions are 


ii.. 

l 


^2 

— = 0 

l 


then 


6^4 =VS- 


2S - 1 


4^(S - 1)tt 2 

This expression also results in a family of curves with * 2 as parameter. 


Center of Gravity at Outboard Quarter Point (Case 4) 

The relations among l and l for this case are 

h..' 

I 4 

> 


The solution to the frequency equation is 


13S 


6 * 3*4 =Vs ' 


4 VS(S - 1)* 2 


(B16) 


(B17) 


(B18) 


(B19) 
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